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Geometric scaling in high-energy QCD at nonzero momentum transfer 
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We show how one can obtain geometric scaling properties from the Balitsky-Kovchegov (BK) 
equation. We start by explaining how, this property arises for the fe-independent BK equation. 
We show that it is possible to extend this model to the full BK equation including momentum 
transfer. The saturation scale behaves like max(g, Qt) where q is the momentum transfer and 
Qt a typical scale of the target. 



Geometric scaling^ expresses the fact that the DIS HERA data depend of the virtuality 
Q 2 and the rapidity Y = log(l/x) only through the ratio Q/Q S (Y) where Q s is an energy- 
dependent scale. On another side, a lot of research have been done during the last ten years 
to obtain a description of the saturation effects in the high-energy limit of QCD, leading to the 
Balitsky-Kovchegov equation^. This equation supplements the BFKlH equation in including 
the contributions from fan diagrams. It has recently been shown in the framework of the impact- 
parameter independent BK equation that one can infer geometric scaling from the asymptotic 
solutions of the evolution equation. 

We shall explain how we have extended these properties of geometric scaling to the case of 
the full BK equation. Our analysis is performed in momentum space where we also introduce 
a new form of the BK equation. We shall discuss the properties of the saturation scale and its 
momentum-transfer dependence. 

Let us thus start our analysis by considering the case of the impact-parameter independent 
BK equation: jV(x, y) = jV(|x — y|). If we go to momentum space using 



l 

2^ 



err 



,ik.r 



M{r) 



dr 



J {kr)M{r) 



o 



J On leave from the fundamental theoretical physics group of the University of Liege. 



0.01 1 100 10000 le+06 le+08 

Rp 



0.01 1 100 10000 le+06 le+08 

Rp 



(a) \og(Rq) = -5 (b) log(Rq) = 2 

Figure 1: Rapidity evolution of the dipole density as a function of p for different values of q (see text). For each 
plot, we show the amplitude for Y varying between and 25 by steps of 2.5. 

the BK equation takes the formal 

d Y M(k) = a X (-d L )M(k) - aM 2 (k), 

where x{l) = 2i/>(l) — VKt) — V^l ~~ l) ls the BFKL kernel and L = \og{k 2 /k^). If we expand 
this kernel to second order around 7 = 1/2, this equation becomes^ equivalent to the Fisher- 
Kolmogorov-Petrovsky-Piscounov equation"! This equation, well-known is statistical physics, 
admits traveling- wave solutions at large rapidities i.e. the solution depends only on Y — v c log(A;), 
where v c is a constant. More generally, if the evolution equation satisfies the following conditions: 

• N = is an unstable fix point and the equation contains a linear term, leading to a growth 
of the amplitude, and a non-linear damping term; 

• the initial condition is steep enough, which is ensured by colour transparency, 

• the linearised equation admits a superposition of waves as solution 

[■c+ioo j 

M M (k,Y) = J ^I ao (7)exp [u(j)Y - 7 log(k 2 /k 2 )] , (1) 
then the asymptotic solution of the full equation takes the form 
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Q 2 s(Y) 



for k 3> ko where the saturation scale Q 2 (Y) grows as fco e VcY . The critical exponent 7c and the 
critical speed v c are given by v c = uj{^ c )/^ c = 0/(7,-), meaning that the group velocity is equal 
to the phase velocity. 

The aim of our studies^^is to show that the properties of having traveling waves as asymp- 
totic solutions holds for the complete BK equation 

cVAA(x, y) = - / d 2 f (x ~ y)2 [AT(x, z) + AT(z, y) - AA(x, y) - AT(x, z)AA(z, y)] . 
7r J (x-z)^(z-yj^ 




It directly appears that, already at the level of the linear terms, this equation is non-local in 
impact parameter. Our first step is thus to move to momentum space 



A/Xk,q) = * / d 2 xd 2 y^e^y^% 
(2ttY J (x - y) z 



After some small manipulations, the equation takes the form 
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d k' M(k, k') AA(k - k', q - k'). 



(3) 



Due to the fact that in this equation the effect of the linear kernel remains local in q, to 
the contrary of what happens in coordinate space, the properties of geometric scaling and the 
formation of traveling waves is better studied in momentum space. 

We now need to check that the three conditions introduced in the previous sections are indeed 
fulfilled in the case of the complete BK equation. The only technical point is obviously to check 
that the solution of the linear equation can be written under the form of a superposition of 
waves. Actually, one can show that the solutions of equation Q, reduced to the linear terms, 
can be put into the form 
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where / 7 (k, q) is known analytically and </>o(7>q) contains all reference to the initial condition. 
Although / 7 (k, q) is a rather complicated function, in the limit \k\ ^> \q\ it simply becomes 
proportional to (k/q)" 2 ^. This means that in the latter limit, the condition (^Q) is satisfied 
with ko = q. Hence, when rapidity increases, traveling waves get formed with the same critical 
properties as in the 6-independent case. The saturation scale is 
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and the tail of the wavefront is then given by 



Af(k,q) ~ 
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Due to the presence of the 0o(7)Q) i n ©> when q becomes smaller than the characteristic 
size Qt of the target, one gets an extra-factor |Qt/<?| 27 - As a consequence, Qt substitutes to q 
as the reference scale fco- The saturation scale Q s becomes is therefore proportional to q or Qt 
whether q is larger or smaller than Qt- 

We have performed^ numerical simulations of the full BK equation in momentum space (J3J. 
For the purposes of these studies, we have introduced R = Qt 1 , the typical size of the target 
and p = k — q/2, canonically conjugated to the dipole size r. The strong coupling constant 
has been fixed to a = 0.2. As clearly shown in figure we have formation of a traveling wave 
when k S> q and at large rapidities. To observe this fact more precisely, we have constructed 
the saturation scale Q s and studied its rapidity evolution. The result is presented in figure [21 
and shows, as expected, goes to v c = 4.886a = 0.9772. 

Finally, figure [2] also shows the q dependence of the saturation scale at large rapidity (Y = 
25). The curve shows that Q s goes to a constant value when q <C Qt and then takes a linear 
behaviour. The very-large-g part of this curve deviates from this behaviour due to the fact that 
the asymptotic regime is not yet reached. 

To summarize, we have been able to obtain traveling waves solutions for the full BK equation. 
For locality reasons, our analysis is done in momentum space. We predict geometric scaling in 
terms of k/(q£l s (Y)) or fc/(Qj'f2 s (y)) whether q is larger or smaller than Qt- The rapidity 
dependence of the saturation scale remains the same as in the 6-independent situation. 

To conclude, numerical simulations would allow for testing these observations phenomeno- 
logically. On theoretical grounds, one can include additional effects such as the running of the 
coupling, the odderon effects or the fluctuation term. 
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